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Abstract. Let P C R'' be a d-dimensional ?i-point set. A Tuerberg-partition of P is a partition 
of P into r sets Pi, . . ., Pr such that the convex hulls conv(Pi), . . ., conv(Pr) have non-empty 
intersection. A point in n[=i conv(Pi) is called a Tverberg point of depth r for P. A classic result 
by Tverberg implies that there always exists a Tverberg partition of size \n/{d + 1)], but it is 
not known how to find such a partition in polynomial time. Therefore, approximate solutions 
are of interest. 

We describe a deterministic algorithm that finds a Tverberg partition of size n/4(d-|- 1)^ in time 
^O(iogd)^ This means that for every fixed dimension we can compute an approximate Tverberg 
^ ■ point (and hence also an approximate centerpoint) in linear time. Our algorithm is obtained by 

combining a novel lifting approach with a recent result by Miller and Sheehy [9]. 
Keywords, discrete geometry, Tverberg theorem, centerpoint, approximation, high dimension 



^ . 1 Introduction 

Let P C M'^ be a d-dimensional point set with n points. In many applications (such as 

j~^ . statistical analysis or mesh generation) we would like to have a way to generalize the one- 

^D I dimensional notion of a median to the high-dimensional point set P. A very natural means 

to accomplish this are centerpoints: a point c € M*^ is a centerpoint for P if every halfspace 

that contains c meets P in at least n/{d + 1) points. Equivalently, for every direction v the 

projection of c onto v splits the projection of P onto v into pieces with no more than dn/{d+l) 

^\ ' points [3]. A classic result in discrete geometry, the centerpoint theorem, shows that there 

j^ ■ exists a centerpoint for every point set [5, 10]. 

However, if we actually would like to compute a centerpoint for a given point set, the 
situation becomes more involved. Helly's theorem implies that the set of all centerpoints is 
given by the intersection of 0{n ) halfspaces [6], so we can find a centerpoint in 0{n ) time 
through linear programming. Chan [1] shows how to improve this running time to 0{n'^~^) 
steps in expectation. He actually solves the harder problem of finding a point with maximum 
Tukey depth. The Tukey depth of a point c' is defined as the minimum number of points in P 
that are met by any halfspace containing c'. A centerpoint has Tukey depth n/{d + 1). If the 
dimension is not fixed, a result by Teng shows that it is co-NP-hard to check whether a given 
point is a centerpoint [11]. In two dimensions, a center point can be found in linear time [7]. 
Since a running time of 0{n'^~^) is not feasible for large d, it makes sense to look for faster 
approximate solutions. A classic approach uses e-approximations [2] : in order to obtain a point 
of Tukey depth n(l/(d + 1) — e) take a random sample ^ C P of size 0{{d/e'^) log{d/e)) and 
compute a centerpoint for A, using the linear-programming method. This gives the desired 
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approximation with constant probability, and the resulting running time after the sampling 
step is constant. What more could we possibly wish for? For one, the algorithm is Monte- 
Carlo: with a certain probability, the reported point fails to be a centerpoint, and we know 
of no fast algorithm to check its validity. This problem can be solved by constructing the e- 
approximation deterministically [2], at the expense of a more complicated algorithm. Nonethe- 
less, in either case the resulting running time, though constant, still grows exponentially with 
d, an undesirable feature for large dimensions. 

This situation motivated Clarkson et al [3] to look for more efficient randomized algorithms 
for approximate centerpoints. They give a simple probabilistic algorithm that computes a 
point of Tukey depth 0{n/{d + Vf) in time 0{(f{d\ogn + log(l/(5))'°s('^+2)), where 5 is the 
error probability. They also describe a more sophisticated algorithm that finds such a point in 
time polynomial in n, d, and \og{l/5). Both algorithms are based on a repeated algorithmic 
application of Radon's theorem (see below). Unfortunately, there remains a probability of 5 
that the result is not correct, and we do not know how to detect a failure efficiently. 

More than ten years later. Miller and Sheehy [9] launched a new attack at the problem. 
Their goal is to develop a deterministic algorithm for approximating centerpoints whose run- 
ning time is subexponential in the dimension. For this, they use a different proof of the center- 
point theorem that is based on a result by Tverberg: any d-dimensional n-point set can be par- 
titioned into r = \n/[d+iy\ sets Pi, . . ., P^ such that the convex hulls conv(Pi), . . ., conv(Pr) 
intersect. Such a partition is called a Twer&erg-partition of P. Any point in 01=1 conv(Pj) must 
be a center point. Miller and Sheehy describe how to find r' = \n/2{d + 1)^] disjoint subsets 
Qi, . . ., Qy.1 of P and a point c' G M'^, such that each Qi contains d+ 1 points and such that 
c S corw[Qi). Hence, c' constitutes an approximate centerpoint for P, and the Qi provide a 
certificate for this fact. The algorithm is deterministic and runs in time vPy^'^^'^i . At the same 
time, it is the first algorithm that also finds an approximate Tverberg partition of P. The 
running time is subexponential in d, but it is still the case that n depends exponentially on 
logd. 

In this paper, we show that the running time for finding approximate Tverberg partitions 
(and hence approximate centerpoints) can be improved. In particular, we show how to find a 
Tverberg partition for P that contains n/4(d + 1)"^ sets in deterministic time d^^^°^'^'n. This 
is linear in n for any fixed dimension, and the dependence on d is only quasipolynomial. 

1.1 Some discrete geometry 

We begin by recalling some basic facts and definitions from discrete geometry [8]. A classic 
fact about convexity is Radon's theorem. 

Theorem 1.1 (Radon's theorem). For every set P C W^ with d + 2 points there exists a 
partition (PijPj) of P such that conv{Pi) D conv{P2) ^ 0. 

As mentioned above, Tverberg [12] generalized this theorem for larger point sets. 

Theorem 1.2 (Tverberg's theorem). Every set P C M'^ with n = {r — l){d + l) + 1 points 
can be partitioned into r sets Pi, P2, . . ., Pr such that n[=i conv{Pi) ^ 0. 

Let P be a set of n points in W^. We say that x € M*^ has Tverberg depth r (with respect 
to P) if there is a partition of P into sets Pi, . . . , P^ such that x G 01=1 conv(Pi). Tverberg's 
theorem thus states that, for any set P in W^, there is a point of Tverberg depth at least 
[(n — l)/[d + 1) + Ij = \n/{d + 1)]. Note that every point with Tverberg depth r also has 



Tukey depth r. Thus, from now on we will use the term depth as a shorthand for Tverberg 
depth. As remarked above, Tverberg's theorem immediately implies the famous centerpoint 
theorem: 

Theorem 1.3 (Centerpoint theorem). For any set P of n points in M there is a point c 
such that all half-spaces containing c contain at least \n/{d-\- 1)] points from P. 

Finally, another classic theorem will be useful for us. 

Theorem 1.4 ( Carat heodory's Theorem). Let P be a set of n points in W^ and x € 
conv{P). Then there is a set of d + 1 points P' (^ P such that x € conv{P'). 

This means that, in order to describe a Tverberg partition of depth r, we only need 
r • (d + 1) points from P. This observation is also used by Miller and Sheehy [9]. They also 
observe that replacing d + 2 by d + 1 points can be done in 0{d'^) time by using Gaussian 
elimination. We denote the process of replacing larger sets by sets of size d + 1 as pruning. 

1.2 Our contribution 

We describe our results in more detail. In Section 2, we present a simple lifting argument 
which leads to an easy Tverberg approximation algorithm. 

Theorem 1.5. Let P be a set of n points in M'^ in general position. One can compute a 
Tverberg point of depth n/2 for P and the corresponding partition in time d^^'n. 

By using this algorithm and combining it with Miller and Sheehy's algorithm, we find a 
point of even higher depth in linear time for any fixed dimension: 

Theorem 1.6. Let P be a set of n points in W^ . Then one can compute a Tverberg point of 
depth n/4(d + 1)^ and a corresponding partition in time 2^^ ^^ ' + d^^' ■ n. 

To the best of our knowledge, this is the first algorithm that computes a Tverberg point 
of polynomial approximation factor in 0{n) time for every fixed dimension. We then improve 
our algorithm to yield an algorithm whose running time is quasipolynomial in d: 

Theorem 1.7. Let P be a set of n points in W^. Then one can compute a Tverberg point of 
depth n/4((i + 1)^ and a corresponding pruned partition in time d ( S") . n. 

2 A simple fixed parameter algorithm 

Next, we present a simple algorithm that runs in linear time for any fixed dimension. For this, 
we show how to compute a Tverberg point by recursion on the dimension. As a byproduct, 
we obtain a quick proof of a weaker version of Tverberg's theorem. 

2.1 The lifting argument and a simple algorithm 

Let P be a d-dimensional point set. A natural way to compute a Tverberg point for P is to 
first project P to some lower-dimensional space, then to recursively compute a good Tverberg 
point for this projection, and finally to try to lift this solution to the higher-dimensional 
space. Surprisingly, we are not aware of any argument along these lines having appeared in 
the literature so far. 



In what follows, we will describe how to lift a lower-dimensional Tverberg point into some 
higher dimension. Unfortunately, this process will come at the cost of a decreased depth for 
the lifted Tverberg point. For clarity of presentation, we first explain the lifting lemma in its 
simplest form. In Section 3.1, we then state the lemma in its full generality. 

Lemma 2.1. Let P be a set of n points in W^. Furthermore, let h be a hyperplane in W^ such 
that the projections of the points in P onto h are pairwise distinct. Let c' € h be a Tverberg 
point of depth r for the projection of P onto h, and suppose we know a corresponding Tverberg 
partition. Then we can find a Tverberg point c ^M. of depth [r/2j for P and a corresponding 
Tverberg partition in time 0{n). 

Proof. For any point p G P, let pr(p) denote the projection of p onto h, and for any Q Q P, 
let pr((5) be the projections of all the points in Q. Let Pi, P2, . . ., Pr '^ P such that pr(Pi), 
pr(P2)) • • •) pr(-fr) is a Tverberg partition for pr(P) with Tverberg point c'. Let I be the line 
orthogonal to h that passes through c'. 

Since our assumption implies c' € conv(pr(Pi)) for z = 1, . . . , r, it follows that i intersects 
each conv(Pj) at some point Xi. Let Qi = {xi-^,Xi^}, i = 1, . . . , [r/2j, be a Tverberg partition 
of xi, . . . ,Xr. Since the points Xi lie on the line i, such a Tverberg partition exists and can 
be computed in time 0{r) by finding the median according to the order along i [4]. Let c be 
some point in the intersection Hil:! conv(Qj).^ We claim that c is a Tverberg point for P of 
depth \r/2\. Indeed, we have 

c G conv(Qj) = convdxjj, Xjj}) C conv(Pj^ U Pi^), 

for 1 < i < [r/2\ . Thus, if we set Qi := Pj^ U Pjj, then Qi, . . ., Q[.r/2\ is a Tverberg partition 
for the point c. The total time to find c and the Qi is 0{n), as claimed. See Figure 1 for a 
two-dimensional illustration of the lifting argument. D 



Fig. 1. Illustrating the lifting lemma in the plane: we project the point set P the line h and find a Tverberg 
partition and a Tverberg point c for the projection. Then we construct the line £ through c' that is perpendic- 
ular to h, and we take the intersection with the lifted hulls of the Tverberg partition (dashed lines) . We find 
a Tverberg partition and Tverberg point c for the intersections. Finally, we group the dashed lines according 
to this Tverberg partition. The resulting convex hulls intersect in the red area, and each of them contains c. 



The proof of Theorem 1.5 is now a direct consequence of Lemma 2.1. 



^ Note that c is either the median or the xt's (if r is odd) or a point in the interval between the median and 
its successor (if r is even) . 



Proof (of Theorem 1.5). If d = 1, we can immediately find a Tverberg point and a corre- 
sponding partition by finding the median c of P [4] and pairing each point to the left of the 
median with exactly one point to the right of the median. 

If d > 1, we project the points onto the hyperplane h : x^ = 0. By general position, this 
results in an n-point set P' C M''"^. We recursively find a Tverberg point and a Tverberg 
partition of depth n/^ for P', and then apply Lemma 2.1 to the result. If desired, we can 
apply the pruning step described in the introduction to the resulting partition, in order to 
ensure that the size of the resulting sets is d + 1. 

The running time for each level of the recursion is dP^^'n, and there are d levels, which 
implies the result. D 

In particular, we obtain a weak version of Tverberg's theorem with a very quick proof. 

Corollary 2.1 (Weak Tverberg theorem). Let P he a set of n points in M . Then P can 
he partitioned into n/2'^ sets Pi, ... , Pn/2^ such that 

n/2'^ 

D 

We now present an alternative, somewhat dual, formulation of Lemma 2.1. For this we 
need some more notation. Let P <Z M. . Recall that a fc-dimensional flat F C M" (often 
abbreviated as k-fiat) is defined as a fc-dimensional affine subspace of M (or, equivalently, 
as the affine hull of k + 1 affinely independent points in W^). A A;-dimensional flat F C M^ is 
called a Tverherg k-flat of depth r for P, if there is a partition of P into sets Pi, . . . ,Pr such 
that conv(Pj) n F 7^ for all i = 1, . . . , r. Note that a Tverberg point in the usual sense is a 
Tverberg 0-flat. 

In the following, we show how to compute a Tverberg (k — l)-flat of depth r/2 from 
a Tverberg fc-flat of depth r for A; > 1. As M'^ is Tverberg d-flat, this again guarantees a 
Tverberg point of depth n/2 . 

Lemma 2.2. Let P he a set ofn points in M . Suppose we have a Tverherg k-flat F of depth r 
for P, as well as a corresponding pruned Tverherg partition. Then, we can compute a Tverherg 
{k — l)-fiat of depth r/2 for P and the corresponding pruned partition in time d^^'n. 

Proof. Let P = l+|J^^ Pi be a partition for the Tverberg k-Hat F. For each i, let Qi be a point 
in the intersection of conv(Pj) with F, and let Q := {qi \ 1 < i < r}. Note that Q can be 
found in d^^^'n time, since we assumed that our Tverberg partition is pruned. Furthermore, 
let hp be a (A; — l)-flat in F that bisects Q {hp can be found by taking the median c' of 
the projection of Q onto a 1-dimensional subflat F' of F and then taking the orthogonal 
complement of F' through c inside F. Let 

Q+ ■.= {q\qeQ,hF-q>Q} = {q^, ..., (?+^ J, 

and Q- ■.= Q\Q+ = {qj^ ,..., q^^^^}. Finally, let Pi := Pi^ U Pj, . 

We claim that hp is a Tverberg (k — l)-flat of depth r/2, with Pi, I = 1, . . . ,r/2 as a 
corresponding Tverberg partition. Indeed, as qf and qj lie on opposite sides of hp, we have 

/ conv(q'+, q~) n /i^ C comr{Pi^,Pj^) = conv(P/), 

for I = 1, . . . , r/2. Since the Pi can be found and pruned in time d^^^'u, the result follows. D 



2.2 Improving the approximation factor 

In order to improve the approximation factor, we will use an easy lemma to bootstrap the 
Tverberg depth. 

Lemma 2.3. Suppose for any m-point set Q (IMr we can compute a point of Tverberg depth 
m/p and a corresponding Tverberg partition in time f{m,d). Let PCM" with \P\ = n, and 
let c € [2,n/p] be a constant, such that n > cp. Then we can find (c — l)p/{d + 1) many 
Tverberg points of depth n/cp for P and their corresponding partitions. These partitions have 
the additional property that each point of P appears in at most one partition. This takes total 
time 

0{j^-{p{n,d) + f{n,d)) 
where p{n, d) is the time for the pruning phase. 

Proof. Let Pi := P. We take an arbitrary subset P[ C P^ with n/c points and find a Tverberg 
point ci of depth n/cp and a corresponding Tverberg partition V'l for P[. Then we prune V'l to 
get a Tverberg partition Vi. Note that this takes time p(n, d)+f{n, d) and that Qi := UzePi ^ 
contains at most n{d + l)/cp points. Set P2 := Pi \ Q2 and continue. 

We repeat this process for a total of a := (c— l)p/{d + 1) iterations, obtaining Tverberg 
points ci, C2, . . . , Cq, with corresponding Tverberg partitions 'Pi, 7^2) • • • 1 "Po- Each partition Vi 
partitions a set Qi ^ P such that the Qi are pairwise disjoint. Furthermore, note that 

\P\> {l-'^^^±^]n> (l-^^^±^]n = ^. 



cp / \ cp 

for z = 1, . . . , a. Thus each q has depth at least n/cp, as desired. D 

For example, by Theorem 1.5 we can find a point of depth n/2'^ and a corresponding 
pruned partition in time d^^^'n. Thus, by applying Lemma 2.3 with c = 2, p = 2'^, we can 
also find 2 /{d + 1) points of depth n/2 ~^^ in linear time, assuming n > 2 ~^^. 

We will also need a lemma that states that sometimes we can multiply the Tverberg depth 
of points. This generalizes a similar lemma by Miller and Sheehy [9, Lemma 4.1]. 

Lemma 2.4. Let P be a set of n points in M , and let P = l+)j^i Pj be a partition of P. 
Furthermore, suppose that for each Pi we have a Tverberg point Ci € M.'^ of depth r, together 
with a corresponding pruned Tverberg partition Vi. Let C := {ci \ 1 < i < k} and c be a 
point of depth r' for C , with corresponding pruned Tverberg partition C. Then c is a point of 
depth r -r' for P. Furthermore, we can find a corresponding pruned Tverberg partition in time 
dOWn. 

Proof. For i = 1, . . . ,k, write Vi = {Qn, ■ ■ ■ , Qir}, and write C = {Di, . . . , Dr'}. For a = 
1, . . . , r', we define r sets Zab as 

^ab '■= {J Qib, 

for 6 = 1, . . . , r. We claim that the set Z := {Zab \ o, = 1, . . . ,r';b = 1, . . . , r} is a Tverberg 
partition of depth r • r' for P with Tverberg point c. Clearly, by definition Z is a partition 



with the appropriate number of elements. It only remains to check that c G conv(Zab) for 
each Zab- We have 

c G conv{Da) C convf (J QiJ = conv(Za6), 

for a = 1 . . .r',b = 1 . . . r. Since certainly \Z\ < n, and since each Zab can be pruned in d^^^^ 
time, the lemma follows. D 

Combining Lemmas 2.3 and 2.4 with the algorithm by Miller and Sheehy, we can now 
prove Theorem 1.6. 

Proof (of Theorem 1.6). If n < 2'^'^^ , we can immediately apply the algorithm by Miller and 
Sheehy to obtain the result. 

Otherwise, we apply Lemma 2.3 as in the remark following that lemma to obtain a set 
C of 2'^ /{d + 1) Tverberg points for P of depth n/2'^~^^ with corresponding pruned partitions 
in linear time. We then use the algorithm of Miller and Sheehy to get a Tverberg point for 
C with depth \C\/2{d + 1)^ with a corresponding partition, in time |Cp('°g'^) = 2'^('^i°S'^). 
Finally, we apply Lemma 2.4 to obtain a Tverberg point and corresponding partition in time 
d^^^'n. The total running time is 2^^'^^°^'^' + d'^^^'n, and the resulting depth is 

n \C\ n 2"^ n 

2^+1 ' 2{d + lf " 2^+1 ■ 2{d + lf " 4{d + lf' 

as desired. D 

3 An improved algorithm 

Finally, we show how to improve our approach through an improved recursion and how to 
obtain an algorithm with running time d^ ^^ ' ■ n. First, however, we describe the more 
general version of Lemma 2.1 we promised above. 

3.1 A more general version of the lifting argument 

We present a more general version of the lifting argument, unifying the previous Lemmas 2.1 
and 2.2. 

Lemma 3.1. Let P be a set of n points in M , and let h Q M. be a k-fiat, such that the 
projection pr{P) of P onto h has n distinct points. Suppose we have a Tverberg point c of 
depth r for pr{P), as well as a corresponding Tverberg partition. Let /i^ be the {d — k) flat 
orthogonal to h that passes through c. Then hjr is a Tverberg [d — k)-flat for P of depth r, 
with the same Tverberg partition. 

Proof. Let pr(Pi), . . ., pr(Pr) be the Tverberg partition for the projection pr{P). It suffices 
to show that conv(Pj) intersects h^ for i = 1, . . . ,r. Indeed, for Pi = {pii,Pi2, ■ ■ ■ ,Piii} 
let c = X]/=i Ajpr(pjj) be a convex combination that witnesses c € conv(pr(Pj)). We now 
write each pij = pr(pjj) + pr^(pjj), where pr-'-(-) denotes the projection onto the orthogonal 
complement h of h. Then, 

as claimed. D 



First of all, this shows how a good algorithm for any fixed dimension improves the general 
case: 

Corollary 3.1. Let 5 > 1 be a fixed integer. Suppose we have an algorithm A with the 
following property: for every point set Q C M. , the algorithm A constructs a Tverberg point 
of depth \Q\/p for Q as well as a corresponding pruned Tverberg partition in time f{\Q\)- 

Then there exists an algorithm such that for any d > 6 and for any n-point set P C R"^ 
we can find a Tverberg point of depth nj p'^' and a corresponding pruned partition in time 
r(i/(51/(n)+dOWn. 

Proof We use induction on k := [d/(5] to show that such an algorithm exists with running 
time k{f{n) + d^^'n). If /c = 1, we can just use algorithm A and there is nothing to show. 

Now suppose k > 1. Let /i C R'^ be a 5-flat in M'^, and let pr(P) be the projection of 
P onto h. We use algorithm A to find a Tverberg point c of depth n/p for pr(P) as well 
as a corresponding pruned partition pr(Pi), pt{P2), . . ., pr(P„/p). This takes time /(n). By 
Lemma 3.1, the {d — 6)-G.at h-^ is a Tverberg flat of depth n/p for P, with corresponding 
Tverberg partition Pi, P2, . . ., Pn/p- For each i, we can find a point qi in conv(Pj) n /i^ in 
time d*^(^). 

Now consider the point set Q = {qi,q2, ■ ■ ■ , Qn/p} ^ ^^- The set Q is a (d — (5)-dimensional 
point set, and since \{d — S)/5~\ = k — 1, we can recursively find a Tverberg point c' for 
Q of depth \Q\/p'^' = nj p'^l and a corresponding pruned Tverberg partition Q in time 
(A; — l)(/(n) + dP^^'n). Now, d is a Tverberg point of depth nj p ' for P, and a corre- 
sponding Tverberg partition is obtained by replacing each point qi in the partition Q by the 
corresponding subset Pj. The resulting partition can be pruned in time dP^^'n. 

Thus, the total running time is {k - \\f{n) + d^^^^n) + /(n) + d^'^^'^n = k{f{n) + d^^^^n), 
and since k = 0(d), the claim follows. D 

3.2 An improved algorithm 

We will now show how to combine the above techniques for an algorithm with a better running 
time. The idea of the new algorithm is follows: using Corollary 3.1, we reduce solving a d- 
dimensional instance to solving two instances of dimension d/2. This can be done recursively. 
Unfortunately, applying Corollary 3.1 reduces the depth of the partition. To fix this, we apply 
Lemmas 2.3, 2.4 and the Miller-Sheehy algorithm to increase the depth again. Details follow. 

Proof (of Theorem 1.7). We prove the theorem by induction on d. As stated before, for d = 1 
the claim is immediate, as in this case the problem reduces to a median computation. 

Thus, suppose that d > 2. By induction, for any ((i/2)-dimensional point set Q CI R'^/^ 
there exists an algorithm that returns a Tverberg point of depth \Q\/A{d/2 + 1)^ and a cor- 
responding pruned Tverberg partition in time d" ^s' "m, for some sufficiently large constant 
a > 0. 

Thus, by Corollary 3.1 (with 5 = d/2), there exists an algorithm that finds a Tver- 
berg point for P of depth n/16{d/2 + 1)^ and a corresponding Tverberg partition in time 

Now we apply Lemma 2.3. The lemma shows that we can find 16{d/2 + 1)^ /{d+ 1) points 
of depth n/32{d/2 + l)^ and corresponding disjoint pruned partitions in time d°^ ^^^ ''^''~^^^' n. 

Let C be the set of these Tverberg points. Applying the Miller-Sheehy algorithm, we can 
find a Tverberg point for C of depth \C\/2{d + 1)^ and a corresponding pruned Tverberg 



partition in time |C'p'-^°^'^''. Now, Lemma 2.4 shows that in additional dP^^'n time, we obtain 
a Tverberg point and a corresponding Tverberg partition for P of size 

n 16(d/2 + if _ n 

2 • 16(d/2 + 1)6 ■ 2(tf + l)2.(d + l) ~ 4(d+l)3' 

as desired. 

It remains to analyze the running time. Adding the various terms, we end up with a time 
bound of 

Since \C\ = dP^^\ we get 

r(n, d) < d"l°gW2)+0{i)^ + (iO{logd)^ < d^l^S'^-^/^n + d^'°^\ 

for a large enough and some /? > 0, independent of d. Hence, it follows that for large enough 
a we have T{n,d) < d°^^°^'^n = d'^^^^^'^'n, as claimed. This completes the proof. D 

4 Outlook 

We have presented a very simple algorithm for finding an approximate Tverberg point, which 
runs in linear time for any fixed dimension. Using more sophisticated methods and combing 
our methods with known results, we managed to improve the running time to d^ ^^ ' ■ n, 
while getting within a factor of l/4((i + 1)^ of the optimum. 

Unfortunately, the resulting running time is still quasi-polynomial in d, and we still do not 
know whether there exists a polynomial algorithm (in n and d) for finding an approximate 
Tverberg point. However, we are hopeful that our techniques constitute a further step in this 
direction and that such an algorithm will eventually be discovered. 

Acknowledgments. We would like to thank Nabil Mustafa for suggesting the problem to 
us and for helpful discussions. 
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